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SUMMARY 


During the first 346 days of the LDEF's almost 6 year stay in space, the metal oxide silicon 
detectors of the Interplanetary Dust Experiment (IDE) recorded over 15,000 impacts, most of which 
were separated in time by integer multiples of the LDEF orbital period (called multiple orbit event 
sequences, or MOES). Simple celestial mechanics provides ample reason to expect that a good deal of 
information about the orbits of the impacting debris particles can be extracted from these MOES, and so 
a procedure, based on the work of Greenberg 1 , has been developed and applied to one of these events, 
the so-called “May swarm”. This technique, the “Method of Differential Precession,” allows for the 
determination of the geometrical elements of a particle orbit from the change in the position of the impact 
point with time. The application of this approach to the May swarm gave the following orbital elements 
for the orbit of the particles striking LDEF during this MOES: a = 6746.5 km; 0.0165 < e < 0.025; 
i = 66°.55; Q 0 = 179°.0 ± 0°.2; cq> = 178°. 1 ± 0°.2. 


INTRODUCTION 


For 346 days after the deployment of the LDEF satellite on April 7, 1984, the tape recorder 
belonging to the Interplanetary Dust Experiment (IDE) stored information on over 15,000 impacts made 
by submicron and larger-size particles on its metal oxide silicon (MOS) detectors 2 . These detectors were 
mounted on trays facing in six orthogonal directions — LDEF ram and trailing edge, the poles of the 
LDEF orbit (north and south), and radially inward (towards the Earth) and outward (towards space). 

The 13.1 second time resolution provided by the IDE electronics, combined with the high sensitivity of 
the MOS detectors and large collecting area (~ 1 m 2 ) of the experiment, conclusively showed that the 
small particle environment at the LDEF altitude of 480 km was highly time-variable, with particle fluxes 
spanning over four orders of magnitude 3 . 


1 Greenberg, R., Orbital Interactions: A New Geometrical Formalism, Astronomical Journal., 87, pp. 184-195 
(1982) 

2 "IDE Spatio-Temporal Impact Fluxes and High Time-Resolution Studies of Multi-Impact Events and Long-Lived 
Debris Clouds", J.D. Mulholland, S.F. Singer, J.P. Oliver, J.L. Weinberg, W.J. Cooke, P.C. Kassel, J.J. Wortman, 
N.L. Montague, W.H. Kinard, LDEF - 69 Months in Space: First LDEF Post-Retrieval Symposium . (NASA CP- 
3134), January, 1992, pp. 517-528 

3 See "LDEF Interplanetary Dust Experiment (IDE) Results", J.P. Oliver et al., this volume. 
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It is highly desirable to learn as much as possible about the orbital characteristics of the particles 
which struck the IDE trays. At a minimum, these characteristics can determine whether the particles 
were interplanetary in origin or debris from a satellite or spent rocket. If the particles can be identified as 
debris, then it becomes possible to determine their parent body, which gives a clue as to which objects in 
Earth orbit are major contributors to the orbital debris population. Unfortunately, the IDE data permit the 
unique determination of only the position of the impacted particle (which is the same as that of LDEF at 
the time of impact), whereas an unambiguous determination of the particle’s orbit requires a knowledge 
of both the position and the particle velocity. The IDE sensors were threshold detectors 4 , triggered by 
any particle with sufficient energy to damage the detector dielectric, and so were rough indicators of 
particle energy, not velocity. It is therefore impossible to use the IDE data to obtain an orbit for a single 
impacting particle. This situation improves, however, for an impacting group of particles which have the 
same orbit In this case, the particles will strike multiple IDE trays, permitting a rough determination of 
the direction of the group’s velocity, which, when combined with the position information, yields a family 
of possible candidate orbits for the particles. The situation improves even more if the orbit of the group 
is such that it encounters LDEF multiple times, for then the change in the LDEF position at the encounter 
times can be used to produce a family of possible orbits, which can be further constrained by the velocity 
direction information. 

Fortunately, most of the 15,000 impacts recorded by IDE occurred in such groups, which we 
term events. These events were of two types — the spikes , which were single, isolated events of high 
intensity and the multiple orbit event sequences (MOES), which were series of events with the events 
separated in time by integer multiples of the LDEF orbital period. The spikes are discussed in another 
paper in these proceedings; here we shall concentrate on the multiple orbit event sequences, as they were 
produced by particles with orbital characteristics such that the group had multiple encounters with LDEF. 
Even though the spikes were generally more intense, the MOES could be quite long-lived, some lasting 
for many days. As discussed in the previous paragraph, it is these MOES which can yield the most 
information about the particles’ orbit. 

Figure 1 is a “seismograph” plot of a typical MOES; time increases to the right along the 
horizontal axis, and the intensities of the events are roughly indicated by the extent of the vertical lines. A 
cursory glance reveals two important bits of information about the particle orbits involved in MOES: 

1) the particle orbits are eccentric; if they were circular, the IDE detectors would register the group 
twice each orbit, as a circular orbit would intersect LDEF’s orbit (which is essentially circular) at two 
points, and 

2) the particles must be “smeared out” along the orbit in some ring-like or torus structure. If the 
particles were concentrated in a “clump”, the encounters with LDEF would not occur at integer 
multiples of the LDEF orbital period, unless the period of the particle orbit was the same as that of 
LDEF, a highly unlikely circumstance. 


4 For details on the IDE detectors, see "Long-term Particle Flux Variability Indicated by Comparison of 
Interplanetary Dust Experiment (IDE) Timed Impacts for LDEFs First Year in Orbit with Impact Data for the 
Entire 5.77 Year Orbital Lifetime", C.G. Simon, J.D. Mulholland, WJ. Cooke, JP. Oliver, P.C. Kassel, LDEF - 69 
Months in Space: Second Post-Retrieval Symposium . (NASA CP-3194), April, 1993, pp. 693-704 


362 



So4h ■ • ■ i i I i I l 


Le4h | | i III I l i • • i i 1 • 


Tr4h 


Figure 1: Typical MOES (not the May swarm). Impacts on the south (So4h), ram (Le4h), 
and trailing (Tr4h) surfaces are shown. The tick marks along the top and bottom are spaced 
at intervals of a LDEF orbital period. 

This information is about all that can be determined from a visual inspection of the MOES in the IDE 
data set. Clearly, it is necessary to develop a technique that will extract additional information about the 
particles’ orbit. We have arrived at such a technique, the “Method of Differential Precession”, which 
shall be summarized and applied in the following pages. 


THE METHOD OF DIFFERENTIAL PRECESSION 


Overview 

The goal of the Method of Differential Precession is to obtain the orbital characteristics of the 
particles which struck the IDE detectors during a MOES by an analysis of the time variation of the LDEF 
position over the series of encounters. This analysis makes use of the fact that the non- sphericity of the 
Earth induces the pole of an object’s orbit to precess, resulting in a cyclic change in the position of the 
line of nodes of the orbit (in the case of LDEF, the period of this precession is approximately 53 days). 
The oblateness of the Earth also causes the line of apsides of the orbit to precess, the point of perigee 
advancing if the orbital inclination is low and regressing otherwise. In general, bodies in different orbits 
will have different rates of these precessions, and should two of these orbits intersect, the differences in 
the precession rates will cause the point(s) of intersection to vary with time. If the characteristics of one 
of the intersecting orbits are known, the migration of the point of intersection may be used to determine 
the precession rates and orientation of the unknown orbit, which then may be used to calculate a family of 
candidate orbits. 

This concept is illustrated more clearly in figures 2 and 3, which depict the geometry of the 
situation with regard to LDEF. Following conventional notation, Ql represents the position of the 
ascending node of the LDEF orbit (which is known) and Q p represents the ascending node of the 
unknown orbit of the impacting particles. The position of the perigee of the unknown orbit is represented 
by ax, LDEF’s orbit is essentially circular (e ~ 10" 4 ) and so has no perigee. The inclinations of the two 
orbits are i L and i p , and ul and u p denote the arguments of latitude of the point of intersection, measured 
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counterclockwise along the orbits from the respective ascending nodes. We are using the argument of 
latitude rather than the more conventional true anomaly, v, due to the fact that one of the orbits is 
circular. In the case of an elliptical orbit, the two quantities are related by v = u - ca Figure 2 shows that, 
for any given time, the known quantities Ql and iL determine the location of the unknown orbit’s node, 
Qp, provided that the inclination of the unknown orbit is specified. Similarly, figure 3 shows that Ql, iL 
and ul determine u p if i p is specified. As time progresses, the orbits will precess at different rates, 
resulting in the movement of the point of impact (intersection), with a corresponding change in ul and u p . 


particle 




Figure 2: Differential precession of the lines of nodes 



Figure 3: Precession of the line of apsides 


Assuming that the inclination of the particle orbit is known, our knowledge of the LDEF orbit 
enables us to compute Q p and u p for each impact occurring in a given MOES. The variation in Q p with 
time directly yields the precession rate of the line of nodes of the particle orbit, Q p , which can then be 
used to construct a family of possible candidate orbits by means of the well-known relation 


• -3/i 

Q = J. 


R. 


la( l-e 2 )] 


cosi p . 


( 1 ) 


where R« is the radius of the Earth and h is the second gravitational harmonic. The semi-major axis and 
eccentricity of the particle orbit are denoted by a and e, while n is the mean motion of the particles. The 
family of candidate orbits will have values of a and e specified by equation (1), and can be constrained by 
the simple fact that any candidate orbit must intersect that of LDEF at some point. Information about the 
direction of the particle velocity obtained from the numbers of impacts on the IDE trays during the 
MOES can also be employed to derive the vector intercept of the particles, which further constrains the 
range of allowed orbits. It should be noted that even though this technique can completely determine the 
orientation of the particle orbit (i p , Q p , and to), the lack of velocity information still prohibits a unique 
determination of the orbit’s size and shape. 

Unfortunately, the inclination of the particle orbit is not known, forcing the adoption of an 
iterative scheme in order to achieve a solution. One starts by assuming a reasonable value for the particle 
inclination, which will enable the determination of the u p ’s and Q p ’s at the times of impact, and, 
consequently, Q p . It is also necessary to obtain the rate of the perigee advance of the particle orbit. This 
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can be done by realizing that, at each impact, the position of LDEF must be the same as that of the 
particle, thus 


a{\-e 2 ) 

y — 

L l+ecos(w p -to) 


( 2 ) 


Equation (2) implies that ii p = co , and we see that the time variation of the argument of latitude of the 
impact point, measured along the particle orbit, is equal to the rate of the advance of the perigee. The 
ratio ti p /G) can now be formulated and compared to the theoretical value, which is given by 



to 


-2 cosip 

5 cos 2 i p -1 


(3) 


If the ratios are not equal, then a new i p is calculated according to Newton’s method or some similar 
scheme, and the process repeated until the values agree. In addition to the assumption of no non- 
gravitational forces, this method also requires that all particles striking LDEF during the MOES share the 
same orbit, which is perfectly reasonable in light of the short duration of each event belonging to a 


MOES. 


Methodology 

Based on the above discussion, one may obtain the family of possible particle orbits by proceeding 
as follows: 

1) Obtain the arguments of latitude (ul) of LDEF at all impact times in the MOE. This is a 
simple matter, given the LDEF orbital elements and an orbit propagation code. 

2) Assume an inclination (i p ) for the orbit of the particles. Good starting values are 30 2 , 65 Q , or 
82° as these are representative of most satellite orbits. 

3) The difference between the right ascension of the particle orbit ascending node (Q p ) and that 
of LDEF (D L ) is determined by the argument of latitude of LDEF at the time of impact (ul), the 
inclination of the particle's orbit, and the inclination of LDEF's orbit (iO- The relevant expressions 


AD = Q -a 

p L 

, —sin AD (4) 

u L = tan ; — — — — 

coti p sin i L - cosA12cosj l 

At each impact time, equations (4) may be solved for AD (and hence, Q p ) via Newton's method or 
some other scheme. 


s Greenberg, R., Orbital Interactions: A New Geometrical Formalism , Astronomical Journal., 87, p. 186. (1982) 
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4) The slope of a line fit to the Q p 's and their associated impact times gives the rate of regression 
of the nodal line of the particle orbit, Q p . 

5) Once the AO's have been determined for the impacts, the corresponding arguments of latitude 
for the particle orbit are found from 

cos u p =cosm l cos AO + sinu L sin AO cos i L 

sinu L sin/ L (5) 


6) In the absence of non-gravitational forces, the semi-major axis (a) and the eccentricity (e) of 
the particle orbit remain constant over time. Therefore, equation (2) requires that u p - 0) must also 
be constant, or « p = to . The slope of a line obtained by a linear regression performed on the u p ’s 

and their associated times yields the progression or regression of the line of apsides, o>. 

7) Compute the ratio Q p /co and compare to the theoretical ratio obtained from equation (3), 

which is a function of only the inclination of the particle orbit. If the two are not equal to within a 
specified tolerance, compute a new i p by means of Newton’s method and repeat steps 3 through 6 
until the values agree. 

8) Use the values of i p and Q p to determine a family of possible candidate orbits in (a,e) space by means 

of equation (1). Constrain the range of potential candidates by imposing the requirements that the 
particle orbit must intersect that of LDEF and must not enter the atmosphere (i.e., the perigee must be 
greater than 200 km). 

Application to the May Swarm MOES 

One of the most prominent multiple orbit event sequences observed by IDE began on May 13, 
1984, and so has become known as the “May swarm.” This MOES can be characterized as being of low 
intensity (~3 impacts per orbit) and long duration, lasting for over 20 days (300 LDEF orbits), with 
several hundred impacts recorded on the IDE trays facing in the LDEF ram direction and towards the 
south pole of the orbit, the majority occurring on the south-facing tray. The long duration of this MOES 
made it an especially suitable choice for analysis by the differential precession technique, the only 
drawback being the low intensity of the events. To avoid contamination by the occasional “random” 
impact, the times chosen were those in which the high sensitivity (0.4 micron dielectric thickness) IDE 
detectors on the south tray recorded multiple impacts within the same IDE clock “tick” (13.1 seconds). 
This resulted in a total of 38 points for use in the analysis, spanning a time interval of some 18 days. 

The procedure outlined in the previous section was then applied to these data, with the inclination 
converging to a value of 66°.55 after only a couple of iterations. Table 1 lists the resulting longitudes of 
ascending node and arguments of latitude of the impact points for the particle orbit, along with the times 
of impact (in decimal days from LDEF deploy) and the LDEF arguments of latitude of the impact points. 
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Impact time 

u L (°) 

Q»(°) 

U D (°) 

40.03233 

221.90 

226.6 

339.7 

40.1638 

227.25 

232.0 

337.5 

40.2922 

215.89 

217.9 

342.2 

40.5546 

223.27 

224.7 

339.1 

40.6842 

218.59 

218.4 

341.1 

40.7499 

221.48 

221.3 

339.8 

40.9452 

219.09 

217.2 

340.8 

40.9463 

224.94 

224.0 

338.4 

41.0106 

219.68 

217.4 

340.6 

41.1417 

223.38 

220.8 

339.1 

41.2063 

219.80 

216.2 

340.5 

41.4682 

224.69 

220.1 

338.5 

41.5339 

226.96 

222.3 

337.6 

41.7942 

223.49 

216.5 

339.0 

41.8599 

225.76 

218.7 

338.1 

41.9906 

226.95 

219.2 

337.7 

42.0554 

224.20 

215.5 

338.7 

42.1856 

222.89 

213.1 

339.3 

43.0342 

224.80 

209.6 

338.5 

43.0343 

225.64 

210.5 

338.2 

43.1635 

218.47 

201.3 

341.1 

43.6219 

227.67 

208.9 

337.4 

43.8179 

228.63 

208.7 

337.0 

43.9480 

226.48 

205.3 

337.8 

44.2092 

228.04 

205.3 

337.2 

44.3403 

230.91 

207.8 

336.2 

45.1233 

231.53 

203.2 

336.0 

45.8412 

232.45 

199.4 

335.6 

46.0377 

236.77 

203.2 

334.2 

47.2122 

236.89 

195.4 

334.2 

47.4733 

237.80 

194.7 

333.9 

47.5380 

235.06 

191.0 

334.8 

47.7988 

234.15 

188.1 

335.1 

50.0854 

250.65 

192.3 

330.6 

50.1505 

249.59 

190.6 

330.8 

50.4769 

250.98 

190.0 

330.5 

52.5007 

257.84 

184.4 

329.4 

55.3086 

271.25 

180.9 

328.7 


Table 1: May swarm times of impact (days from LDEF deploy) and the corresponding LDEF arguments of latitude, 
with the values of the longitudes of the ascending node of the particle orbit and the particle arguments of latitude 
for i p = 66°.55. 

The two linear regressions (see figures 4 and 5), involving the impact times, Q p , and u p , yielded 
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Q = 179°.0±0°.2 (Initial longitude of ascending node for the particle orbit) 
pO 

n = -3°.26±0°.05 day' 1 
P 

C0q = 178° . 1 ± 0° .2 (Initial argument of perigee for the particle orbit) 

6) = -0°.85±0°.05 day’ 1 

These four quantities, along with i p , uniquely specify the orientation of the particle orbit at any given time. 
Note that the initial value of the argument of the perigee indicates that these particles are striking LDEF 
near apogee, a somewhat surprising result. 




38 40 42 44 46 48 50 5 2 54 56 58 

Days from LDEF Deploy 


Figure 4: Linear fit to determine nodal line properties 


Figure 5: Linear fit to determine apsidal line properties 


Next, the precession rate of the particle orbit line of nodes was used in equation (1) to determine 
the family of possible candidate orbits. These results are displayed in figure 6. Note that the semi-major 
axis varies little with the eccentricity; in this case, the variation in a is so small that we could confidently 
set a = 6746.5 km, regardless of the eccentricity. The dual requirement that the candidate orbits have 
perigees of greater than 200 km in altitude and intersect the LDEF orbit placed strict limits on the 
allowed values of the eccentricity, which must lie in the range 0.0165 < e < 0.025. 

One of the candidate orbits (e = 0.017) was then chosen for a series of checks on the results of the 
method. The first check involved the computation of the particle velocity of impact over the duration of 
the May swarm. These velocities were then resolved into components along the LDEF body axes in 
order to determine the impact speeds on the IDE trays. For this particular orbit, only the south tray and 
the ram-facing tray were struck, with the south impact speed being larger than that for the other tray (see 
figure 7). This is in good agreement with the IDE observations of the May swarm, in which these same 
two trays recorded large numbers of impacts, with the south tray receiving the most hits. The second 
check consisted of a comparison of the sky track of the points of closest approach between the two orbits 
to the sky positions of the individual impacts comprising the May swarm. As can be seen from figure 8, 
the agreement is excellent, with the sky track of close approach passing neatly through a diffuse band of 
impact positions. 
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Figure 6: Candidate orbits for the May swarm 
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Figure 7: Particle impact speeds along IDE tray normals for test orbit. 

In summary, it would seem that the particles impacting LDEF during the May swarm 
MOES have an orbit that can be characterized by the following parameters: 
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Figure 8: Sky track of close approach between the test orbit and that of LDEF. The particles are moving in a northerly 
direction, whereas LDEF is moving along its orbit from left to right At the onset of the May swarm, the impacts are 
located at the position labeled “Onset”, with the impact positions gradually moving towards the lower right as time 
progresses. 
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Qpo = 179°.0±0°.2 
O p = -3°. 26 ± 0°.05 day’ 1 

cub = 178°.l ±0°.2 
c6 = -0°.85 ± 0°.05 day 1 


ip = 66°.55 
a = 6746.5 km 

0.0165 < e <0.025 


CONCLUSION 


It has been pointed out numerous times in the literature that the surface area to mass ratio 
of a micron-sized particle (the size of many of the LDEF impactors) is large, thus causing the 
particle to experience significant perturbations from forces such as radiation pressure and 
atmospheric drag. Indeed, numerical calculations indicate that radiation pressure can cause a one 
micron diameter particle, initially in an orbit similar to that of LDEF, to enter Earth’s atmosphere 
after only a very few (<10) orbital revolutions. Such calculations leave us hard-pressed to explain 
how a MOES like the May swarm, which we postulate to be caused by a ring of micron and sub- 
micron particles, can persist for many days. The only reasonable explanation is that the ring must 
be replenished by debris from some source during the time spanned by the MOES. 

To obtain the geometrical characteristics of this ring, the technique of differential 
precession looks at the time evolution of the point of intersection with the LDEF orbit. It does 
not matter that the particles exist in the ring for only a short time; the only requirement is that the 
orbit shared by the particles at the times of impact with LDEF be similar. In general, this orbit 
would not be the same as that of the source of the debris particles, for non-gravitational forces 
would have rapidly acted to alter the particles’ orbit from that of the parent body. It should be 
realized that if the parent body (whose orbit is presumably stable) continually produced particles 
of similar properties, these particles would have experienced the same perturbations as their 
predecessors and would therefore have undergone a similar orbital evolution. If any of the future 
orbits intersected that of LDEF, a MOES would have been observed by IDE, this MOES lasting 
as long as the source produced particles, or until the geometry of both orbits changed such that 
there was no longer a point of contact. 

The IDE data set is rich, with many MOES of varying characteristics that await analysis by 
some procedure. The method of differential precession is such a technique, one that appears to be 
able to extract a good deal of information about the particle orbit involved in a MOES. We fully 
expect that its application to the other MOES will not only shed some light on possible sources of 
orbital debris, but will also yield quite a few surprises. 

This research was supported by NASA grant NAG-1-1218. 
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